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Abstract 

The problem of compressive detection of random subspace signals is studied. We consider signals 
modeled as s = Hx where H is an A x iT matrix with K < N and x ~ We say 

that signal s lies in or leans toward a subspace if the largest eigenvalue of HH^ is strictly greater 
than its smallest eigenvalue. We first design a measurement matrix $ = comprising of two 

sub-matrices and where projects the signal to the strongest left-singular vectors, i.e., the 
left-singular vectors corresponding to the largest singular values, of subspace matrix H and projects 
it to the weakest left-singular vectors. We then propose two detectors which work based on the difference 
in energies of the samples measured by the two sub-matrices and and provide theoretical proofs 
for their optimality. Simplified versions of the proposed detectors for the case when the variance of noise 
is known are also provided. Furthermore, we study the performance of the detector when measurements 
are imprecise and show how imprecision can be compensated by employing more measurement devices. 
The problem is then re-formulated for the generalized case when the signal lies in the union of a finite 
number of linear subspaces instead of a single linear subspace. Finally, we study the performance of the 
proposed methods by simulation examples. 
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I. Introduction 

The topic of Compressive Sensing (CS), where the effort is to draw inferences about the sparse 
signals based on compressive samples, has been one of the most attractive topics in the area of signal 
processing for more than a decade. Compressive Sensing states that if a signal is sparse in some known 
bases, then it is possible to reconstruct it from a relatively few linear projections ||T1-||3||. While most 
of the efforts in the area of CS has been devoted to the problem of estimation of a sparse signal from 
possibly noisy compressive measurements, there have been also some efforts in accomplishing other 
signal processing tasks such as Compressive Detection Compressive Classification 121, lfT2 - 

m, etc. The goal of this paper is to study the problem of detecting a random subspace signal from its 
compressive measurements. 

Signal detection ifTTl is an important task in Statistical Signal Processing where the goal is to decide 
on the presence/absence of a signal rather than estimating it. Signal Detection has applications in many 
areas of engineering and science, including wireless communications, radar and sonar, bioinformatics, 
etc. The literature on detection from classical uncompressed samples is very vast and well-studied; see, 
e g-, ifT/l - lT^ and references therein. With the advent of compressive sampling methods, like many 
other fields, there have been additional efforts to tailor the existing signal detection techniques to the 
case where the samples are taken at compressive rates. The first criterion that a signal has to meet to 
qualify for being detectable based on compressive samples is to be sparse in some known basis, or 
in a broader sense to lie in a low-dimensional subspace of the higher-dimensional ambient space. The 
applications of such compressive detection techniques are in the scenarios where uncompressed sampling 
is difficult, expensive, or even harmful, but because of the existence of a structure in the signal we are 
still able to sample the signal in compressive rate and then detect it based on the compressive samples. 
Examples of such scenarios are compressive spectrum sensing of OFDM signals 1201 . iflTl . radar sensing 
II 22 I - II 22 . hyperspectral imaging |[26l . ultrasound imaging ll27l . etc. We remark here that the hardware 
implementation of specific compressive samplers is beyond the scope of this paper and we only focus on 
finding the optimal mathematical model or framework of such samplers, however examples of hardware 
implementations in various applications can be found in 1221 . Il271 - ll33l . 

In this paper we address the problem of compressive detection of signals that lie on or close to a low 
dimensional linear subspace. We start with the simplest case when the signal is drawn from a known 
low-dimensional subspace and propose two tests for detection of signal from compressive measurements. 
The main advantage of the proposed algorithms over the existing methods is that they do not need 
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the knowledge of the noise variance to work. However, we also provide the simplified versions of the 
proposed methods for the case when the noise variance is known. We then study the effect of imprecise 
measurements and show that the imprecision can he compensated hy employing more measurement 
devices. Then the proposed method is re-formulated for the case when the signal is lying on the union 
of a finite number of known linear suhspaces with each of suhspaces having a certain prohahility of 
being the subspace on which the signal truly lies. In other words, in this case, instead of having the 
exact knowledge of the subspace, we know the set of all subspaces from which the signal may emerge 
together with the corresponding probabilities that the signal may lie on each of them. We will introduce 
the optimum sampling strategy in this case and show that although the performance falls by increasing 
the number of possible subspaces, it remains reasonable as long as the cardinality of the set is much 
smaller than the ambient dimension. 

Related work: In 0|-l|6l, ifT^ authors studied the problem of compressive detection, but they formu¬ 
lated the problem for deterministic agnostic signals rather than random subspace signals. Besides, their 
methods did not take into consideration the unknown noise variance scenario. The method in Q designs a 
matched subspace detector for subspace signals. The method assumes that the exact knowledge of signal 
is available for detector design. They extend their work to unknown signals in but like the previous 
ones they need the variance of noise to design the compressive detector. The work of ifTOll provided upper 
bounds for probabilities of false alarm and mis-detection for deterministic signals while the variance of 
noise is again assumed known. Finally, in ifTTl authors studied the problem of compressive detection of 
random subspace signals but similar to the previously mentioned works the noise variance was assumed 
known and the signal was assumed agnostic. 

Notations and Mathematical Preliminaries: Throughout this paper, all quantities are assumed real¬ 
valued while matrices and vectors are denoted by capital and small boldface letters, respectively. = 
denotes the equality and = denotes the definition. Kx is reserved for statistical expectation operator with 
respect to the random variable X, and (8) denotes the Kronecker product. Ip, Ipg and Op^g represent 
P X P identity matrix, and P x Q all-one and all-zero matrices, respectively. The set of real numbers 
is denoted by M, the set of real-valued M x 1 vectors is denoted by the set of real-valued M x N 
matrices is denoted by and for an arbitrary matrix A G [A]jj denotes the (i,j)-th entry 

of the matrix. For a set Ad, cardinality of the set is shown as |A4|. For a matrix A we show its i-th 
eigenvalue by Ai(A). We also use the notations Amax(A) and Amin (A) for denoting its maximum and 
minimum eigenvalues, respectively. For a vector x, its sub-vector containing entries from i to y > i is 
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denoted by x 


Paper Organization: The rest of this paper is organized as follows: first, in Section HIl we explain the 
basic model of a linear subspace signal and we introduce the sampling strategy that we choose for our 
compressive detector. Then based on the introduced sampling strategy, in Section JII] we introduce two 
compressive detectors for detection of the signal based on compressive samples and provide theoretical 
results about their optimality. We also briefly study the simplification when the noise variance is known 
as well as the case of subspace interference. Next, in Section |IVl we study the effect of imprecise 
measurements where, e.g. because of hardware design limitations, the measurement matrix may deviate 
from the ideal one introduced in Section JII] In Section jVj we extend our design to the case where instead 
of knowing the exact linear subspace on which signal lies, we have only more coarse knowledge about 
the signal location in the ambient space in the form of knowing all the possibilities for the true subspace 
together with their corresponding probabilities. The performances of the proposed methods are studied 
through computer simulations in Section |Vl| Finally we conclude the paper in Section I VII I 


II. System Model 

Consider first a noiseless subspace signal of the form 

s[n] = Hx[n], n = l,...,Nb, (1) 

where x[n] ~ Af{OK,i,<^x^K) and H is an x iT full-rank deterministic matrix with K < N. We denote 
the N nonnegative eigenvalues of HH^ by > /?| > ... > > 0 and we further assume that 

Pi > P%- (2) 

Inequality (|2]) implies that the signal energy is not distributed uniformly over all dimensions, or in other 
words, the signal lies in (or at least leans toward) some subspace of the ambient space. 

In a classical uncompressed scenario, by processing a sequence of noisy observations y[n],n = 

1.. .. ,Nb, the signal detection problem refers to the following hypothesis testing 

I ^o:yN=w[n], n = 1,2,... ,iVfe, 

\ T-ii : y[n] = s[n] + w[n], n = 1,2,..., iVfe, 

in which w[n] ~ AA( 07 v,i, Iat) refers to observation noise. However, as mentioned in the previous section, 
obtaining uncompressed observations is sometimes difficult or expensive, and therefore we have to perform 
the signal detection task through some noisy compressive observations z[n] = [zi [n], 22 [n],..., zm [n]]"^,n 

1.2.. .., Nb, rather than unavailable uncompressed samples y[n], n = 1,2,..., Nb. Here Zm[n], m = 
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1, 2,..., M denotes the output of m-th compressive measurement unit. Thus, we redefine the hypothesis 
testing problem as 

I ^o:zN=w'^Ni n = 1,2,..., 

I : z[n] = $s[n] + w'^[re], n = 1, 2,..., A^'b, 
where $ = [0^, (^> 2 ,..., 0^^]^ € is compressive measurement matrix with M N and 

w'^[n] = [<?!)fwi[re], (?!)^W 2 [n],..., 0 ^wmN]^ where ~ A^( 07 v,i, Iat) represents the uncom¬ 

pressed receiver noise at the input of the m-th compressive measurement unit. We have used different 
indices for noise vectors as we assume that these measurement units work independently and posses 
independent noise at their inputs, which means that w^i [n] and [^t-] are statistically independent for 


mi / m2. 

This problem was addressed in lO, lU in presence of noise with known variance and for deterministic 
signals. In this paper, we design compressive detectors without assuming prior noise knowledge. To this 
end, we partition the M measurement devices to two sets or, in other words, we consider the measurement 


matrix $ to consist of two sub-matrices € 


hMixN 


and G 




as 


$ = 




(5) 


where Mi -|- M 2 = M. The idea behind having two sets of samplers, and $o> is to get two sets of 
measurements 

Zs[n] = z[n] , (6) 

and 

z„[nl = z[nl , (7) 

Mi + l-.M 


whose selected statistics are identical under null hypothesis but different under alternative hypothesis. 
The statistics of the two sets of measurements can then be compared to decide whether we should 
accept the null hypothesis Tfo or reject it in favor of alternative hypothesis Tfi. In the next section, we 
provide two designs for and $0 that fulfill this goal and then propose detectors for each design to 
perform the hypothesis testing problem of ([D) based on two sets of compressive measurements z^ [n] and 
Zo[n],n = 1,2,... ,W'b- 


III. Compressive Subspace Detection 

Based on the system model of the previous section, in this section we propose two detectors for 
performing the hypothesis testing problem of ([111. Both detectors are composed of two sets of measure¬ 
ment devices, represented by and $ 0 . In the sequel, we introduce these two detectors, derive their 
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performance measures, and prove their optimality. The numerical study of the performance of the two 
proposed detectors will he provided in Section through simulation experiments. 


A. Maximally-Uncorrelated Compressive Detector 

Assume that the SVD of matrix H G is given hy 

H = UDV'^, (8) 

where U G and V G R^^^ are orthogonal matrices, and D G R^^^ is a diagonal matrix. 

Consider positive integers Mi and M 2 such that M = Mi + M 2 . Let us denote the first Mi columns 
3y Us_and its last M 2 hy Uq. We propose the following design for M x N measurement matrix 

with suh-matrices and $0 defined as 


of U 
$ = 


*0 


and 






1 7, 

_T Tr 


T 

J-O Wq , 


y/M 


(9) 

( 10 ) 
in (O and 


where and To are Mi x Mi and M 2 x M 2 orthogonal matrices, respectively. The term 
([Toll is to guarantee that the received energy at the output of the compressive sampler is independent of 
the number of samplers. This makes the study of the effect of the number of samplers on the performance 
fair as otherwise it is obvious that the bigger is the number of samplers the better is the performance. 
We remark again that the M rows of the measurement matrix $ represent M measurement devices 
equipped to collect samples. Hereafter, we use the terms measurement matrix and measurement devices 
interchangeably. 

The reason for choosing this design for measurement matrix is that the rows of span the Mi- 
dimensional subspace which contains the highest amount of energy of signal s = Hx, or in other words, 
is the Ml-dimensional subspace along which the variance of s is maximized. Furthermore, the rows of 
$0 span the M 2 -dimensional subspace which contains the lowest amount of energy of signal s, or in 
other words, is the M 2 -dimensional subspace along which the variance of s is minimized. The difference 
between the energy of signal taken at the output of these two sub-matrices can then be exploited for 
signal detection. The following theorem, justifies the above discussion. 

Theorem 1. The measurement sub-matrix in (|9l) is the solution to the following optimization problem 


arg max Ex,-v 




S.t. M^s^s =I 


'-Ml, 


( 11 ) 
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In other words, the rows of in (|9l) represent the set of Mi uncorrelated measurement devices that 
maximize the expected value of total energy (or the total variance) in their outputs. 

Correspondingly, the measurement sub-matrix in diOl) is the solution to the following optimization 
problem 

argminEx,w(||zo||i), s.t. = Im^, (12) 

In other words, the rows of in diOl) represent the set of M 2 uncorrelated measurement devices that 
minimize the expected value of total energy (or the total variance) in their outputs. 

Proof: We first notice that 

IEx,w(||zs||i) = Ex(x'^H^$f$,Hx) + ^^ 

= cj2trace(HH^$f$,) + ^^^ (13) 

Ni 

Noticing that the second term in the right-hand side of (fT3l) is independent of measurement design and 
H^H and are both symmetric matrices and that Mi $ J in (fTTI) has Mi eigenvalues equal to 1 

and N — Mi eigenvalues equal to zero, the proof of ® is easily concluded hy employing von Neumann 
trace theorem Il34l P 11.4.5], |[35l Theorem 6.77]. The proof for (fTOl ) is similar. ■ 

Hereafter, we call and $0 as max-energy sampler and min-energy sampler, respectively, and call the 
measurement design in (l9ll and (fTOl) Maximally-Uncorrelated Compressive Detector as it is the optimum 
design when the measurements are uncorrelated. 

Next, based on the above measurement design we introduce a test statistic for carrying out the 
hypothesis testing problem of (01). The following theorem states this test and derives its performance 
measures. 


Theorem 2. Consider the following test for performing the hypothesis testing problem given in 



n=l 


(14) 


where Zs[n] and Zo[n] are as in and dZl), respectively. The Probability of False Alarm, PfA, for this 
test is computed as 


PpA — QF{MiNt,M2Nt)i'y)^ 

and the Probability of Detection, Pd, is bounded as 

Pd, lb < Pd < PD,uh 


(15) 


(16) 
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where PD,ib <ind PD,ub CL^e, respectively, the lower bound and the upper bound for the probability of 
detection which will take the following formulas 


Pd, lb — QF(MiNb,M2Ni)i'nibl), 


(17) 


and 


PD,ub — QF{MiNb,M2Nb)iVub'y), 


(18) 


where pn, = , Pub — QF{di,d 2 )i^) probability of the F-distribution 

with parameters (degrees of freedom) di and d 2 at point x. 


Proof: To prove the theorem, we first notice that Zo[n] has distrihution 

Ar(0, ^Ims) under 7^0, 

Zo[nj ~ < (19) 

[ A7(0, ]^(crglM 2 + o-^ToAoT^)) under Pi, 
where Aq is an M 2 x M 2 diagonal matrix whose (i, f)-th diagonal element is 
On the other hand, the distrihution of z^ is 

A7(0, ^Imi) under Pq, 

ZsN ~ (20) 

A7(0, A(cr2j^^ _P cr2T^AsTf)) under Pi, 
where A* is an Mi x Mi diagonal matrix whose (f,f)-th diagonal element is pf. 

From (I 20 I) we have 

X; Z5[n]^z^[n] 

2 ,.. -~ >under Pq, (21) 

CTqIIVI 

where denotes chi-squared distrihution with d degrees of freedom ifTTI . On the other hand, under Pi 
we have 


Nb 


^Z 5 [n]^A^ ^Zs[n] ~ under Pi. 


( 22 ) 


n=l 


where A^ = + c^T^AsTJ). From Rayleigh-Ritz theorem |[36l Chapter 8] we know that 


-^minCAg )z^ Zg f- Zg A^ Zg f- Amax(Ag )Zg Z 


A-1 




(23) 


where A^in(A, i) = 


and Amax(A 5 ^) — ^ — 


M 


Amin (As) of+Tfp^' 
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Similarly, from (fT^ we have 


Nb 


Zo[n]^Zo[n] 


n=l 


a^M 


under T-Lq, 


and 


Nb 


^Zo[n]^Ao ^Zo[n] ~ under ^i. 


where Aq = ^^(crolMa + u’^ToAoT^). We also have 


n=l 
Sf-p A 'T'T^ 


Amin(Ao )Zo Zo < A^ Zq < Amax(Ao )z^ Z, 

where Ainin(Ao l) = a„J(A„) = -^max(Ao - a-+aipi, ' 

From above computations, under "Hq hypothesis we will have 


^ z^nj^z^n] 

n=l 

"^5 

^ Zo[n]'^Zo[n] 
n=l 


At,MaAb,under "Ho- 


Under "Hi hypothesis we have 


Ab 




E z,[n]^A7iz,[n] 


Ab 


-^MiAb,M2Ab 


yy z<,[n]'^z<,[n] 
('^o+o'xPn)/A^ >^=1 


E Zo[n]^Ao^Zo 


n 


and similarly 


yy;^zs[n]^z4n] 

('^o+o'xPmJ/A^ 


yy z<,[n]'^z<,[n] 

n = l 


> 


Ab 

E z,[n]^A7iz,[n] 

72 = 1 


FMiNb,M2Nb 


E Zo[n]^Ao^Zo 


re 


(o' 0 +'^xPn-M 2 + i)/^ 

From (l27l) . (|2^ . and (l29l) the proof is concluded. 

Remark 1: The Q-function of F-distrihution is expressed as ll37l 

QF(di,(i2)(^) = 1 -^1" -a)) 

d^x + ti2 ^ ^ 

where Ix{a,b) is the regularized incomplete Beta function defined as 

B{x-, a, b) 


Ix{a,b) = 


B{a,b) ’ 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


where B{x\a,b) = Jq ^(1 — t)^ ^dt is the incomplete Beta function and B{a,b) = B{l;a,b) is the 


Beta function. 


□ 


December 31, 2015 


DRAFT 


















10 


Remark 2: From Theorem |2] it is clear that for a given Mi and M 2 the proposed detector is unbiased 
|[T^ Chapter 4] if p^-M^+i < PMi- Furthermore, if p\ > then there is always a design as suggested 
in this section with certain Mi and M 2 which delivers an unbiased detector. □ 

Theorem [T] showed that the proposed design of (|9l) and (fTOb is the one that provides the largest difference 
between the expected energies of signals and Zq under orthogonality condition. Then in Theorem |2] 
a test for detecting random subspace signals based on this design was proposed and the probabilities 
of false alarm and detection for the proposed detector were derived. An important measure in statistical 
hypothesis testing is to show that the proposed test is the most powerful test in the sense that it has 
the highest probability of detection, P^), for a fixed probability of false alarm, PpA- In the following a 
corollary is provided to show that among all orthogonal designs, the one proposed in (|9l) and (ITOl) delivers 
the highest probability of detection for the detector in (fT4b . 


Corollary 1. Among all orthogonal designs satisfying = -^Im with fixed number of measurement 
devices {Mi, M 2 ), the design of ^ and dTO]) provides the highest PD,ib and PD,ubfor a given PpA- 


Proof: Consider an arbitrary orthogonal measurement design as $ = $0 ]^, where G 

^MixN ^ j-Q jf choosc them as in (|9l) and (ITOb then PD,ib 

and PD,ub nre maximized for a given PpA. To this end, first by taking the same steps as in the proof 
of Theorem m it is easy to obtain the Probability of False Alarm as PpA = QF{MiNb,M 2 Nb){l)- Then 
for fixed PpA, Mi, M 2 , and Nh fhe lower and upper bounds on fhe Probabilify of Defection, PD,ib < 


Pd < PD,ub, can be respectively expressed as Pd^w = QF{M^N,,M 2 N,){mbQplM^N,,M 2 NA^^FA)) 


and 


PD,ub — QF{MiNb,M 2 Nb){hubQp(^]if_^]\f,^^ju^]\l^'^{PFA)), where 


mb = ^ 






(32) 


and 

Therefore the optimum design for and $0 is the one that minimizes fin, and 7)^6 in (l32l ) and (|3^ . which 
in turn is the one that minimizes Ai ^oH) and Ams ^oH) and maximizes Ai J^sH) 

and Ami(H^^J^sH). The proof is then easily concluded by employing Poincare Separation Theorem 
041 Theorem P 10.4.2]. ■ 

Remark 3: The proposed detector can be thought of as an F-test detector. F-test 0^ is a statistical 
test for comparing the variances of two normal populations. The proposed detector can be seen as an 
F-test detector where the two populations, i.e., outputs of the two sets of measurement devices, have been 
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designed so as to impose maximum difference to their variances, which optimizes the performance of the 
F-detector since it works based on the difference between variances. This is in fact why the optimizations 
in Theorem [U and Corollary [B which are trying to optimize different cost functions, result in identical 
solutions. □ 

While Theorem |2] proposes a test statistic for signal detection and provides a proof for the performance 
of the proposed detector, it does not provide us with the optimum number of equipped samplers. Because 
of the complicated form of the tail probability of F-distribution and regularized incomplete Beta function, 
finding a closed form solution for Pd as a function of Pfa, Mi, M 2 , and is overly cumbersome. 
However, for large degrees of freedom, say larger than 100, it is possible to provide an analysis based 
on the following Normal approximation of the F-distribution |[3^ Section 12.4.4] 


QF{dud^)ix) ~ Qz{x - jl/a) 


(34) 


where is the tail function of the standard Gaussian distribution, 

d2 


F = 


d2 -2’ 


and 


Then, for our problem, we can write 


a = fii 


l2{di+d2-2) 

di(d2-4) 


(35) 


(36) 


Pd, lb = QzivibdQz ^(^’fa) + riibfj), (37) 

and similarly for PD,ub> where di = aMN^, d 2 = {1 — a)MNf,, and a = Mi/M. From (1^ and (l36l) it 
is clear that for big values of MNi,, both jl and d are decreasing functions of MNi,, and therefore PD,ib 
is an increasing function of both M and iVfe. This can be seen also from Figure [T] which illustrates an 
example of how the performance of the compressive detector changes with the total number of sampling 
devices M and the number of temporal measurements N^. In this example, we set pi = pM^ = 1 and 
Pn = PN-M2+1 = 0 to have Pd = PD,ib = PD,ub- Besides, we choose Mi = M 2 = M/2, Pfa = 0.05, 
and cj^ = (To = 1. The probabilities of false alarm and detection are then calculated from (fTSl) and ([TV]) . 
Notice that from a practical point of view an important factor which limits utilizing high values of Mi 
and M 2 is the hardware cost and complexity as it is proportional to the number of samplers which is in 
turn determined by Mi and M 2 . 

Furthermore, for Ni, ^ 1, one can apply a further approximation on (l35l) and (|3^ to obtain 




1, 


(38) 
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Fig. 1. Probability of Detection as a function of M and Ni,. Here we have Mi = M 2 = M/2. 


and 


CJ 


(39) 


a(l — a)NbM 

Equation (l39l) then shows that the performance is optimized for a = 0.5, or in other words when we 
choose Ml = M 2 . This means that when designing a Maximally-Uncorrelated Compressive Detector, 
the optimum choice is to divide the hardware budget equally between max-energy samplers and 
min-energy samplers $o- This can be seen also from Figure [2] which shows that for a fixed number of 
measurement devices M = 16, the maximum indeed occurs practically in the middle, which means that 
the optimum choice is to allocate half of the measurement devices to $5 and the other half to 


B. Fully-Correlated Compressive Detector 

In this part, we study the effect of increasing the Signal-to-Noise Ratio (SNR) at the expense of 
uncorrelatedness of measurement devices. For a certain measurement device, the effective noise energy 
can be decreased (and therefore the SNR can be increased) by repeating the measurements and averaging 
over all of them. Collecting the multiple measurements one after the other over time for enhancing the 
effective SNR is not however a suitable solution here, firstly because the signal itself is random too, and 
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Fig. 2. Probability of Detection as a function of Mi and Nb for a fixed total number of measurement devices M — 16. 


secondly because it increases the latency which is an issue in many applications such as radar or wireless 
communications due to the relatively quick changes of communication channel characteristics. 

However in our proposed measurement framework, improving the effective SNR can be done by 
employing Mi identical measurement vectors instead of Mi x N measurement matrix in ® and M 2 
identical measurement vectors instead of M 2 x N measurement matrix $0 in (fl^ and then averaging 
over all measurements. This will not increase the latency because the identical measurements are taken 
in parallel. In other words, matrices and $0 in (O and (fTOl) should be respectively redefined as 

1 


® uf, 


and 


y/M 

1 


<S) U^, 


/m 


(40) 


(41) 


where ui and uat are the left singular vectors of H corresponding to the pi and p^, respectively. The 
following theorem is to show the significance of fhe measurement design in (l40l) and (|4T]) . 

Theorem 3. The measurement sub-matrix in d40l) is the solution to the following optimization problem 


arg maxEx.wdIzsE) 

$ 


(42) 
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In other words, the rows of in d40D represent the set of M\ measurement devices that maximize the 
total energy (or the total variance) in their outputs. 

Correspondingly, the measurement sub-matrix (EZ) the solution to the following optimization 

problem 

argminEx,w(||z5||i)- (43) 

$ 

In other words, the rows of (IS]) represent the set of M2 measurement devices that minimize the 

total energy (or the total variance) in their outputs. 

Proof: We first notice that 

IEx,w(||zs||i) = Ex(x'^H^$J$,Hx) + ^^ 

= a2trace(HH^$f$,) + ^^ (44) 

Since the second term on the right-hand side of (l44l) is independent of measurement design, the optimum 
design for is the one which maximizes the first term. Noticing that H^H and are both 

symmetric, from ll34l P 11.4.5] we have 

Ml 

trace(HH^$,$f) < (45) 

i=l 

Unlike Theorem [T] here we have now no constraint on matrix other than its rows having norm IjsfM 
which means that trace($s$J) = Xi{^s^J) = Mi/M. Therefore the design of can he done 
in two steps. In the first step we find the eigenvalues that maximize the right-hand 

side of inequality (l44l) . From BOl 4.B.7 and l.A] the maximum is attained when Ai($s$J) = 1 and 
{Aj($s$J )}^2 = 0- III Ills second step we simply employ |[34l P 11.4.5] to conclude the proof of (l40l) 
being the solution to (l4^ . The proof for (HTl) is similar. ■ 

After designing the samplers as in (l40l) and (|4TI) and collecting measurement vectors [n] and Zq [n ], 

the effective energy of noise can then be reduced by averaging over all elements of z<j[n] and Zo[n] to 
obtain 

= llVl-i,MiZs[n], n = l,...,iVb, (46) 

lVi\ 

and 

^oN = ^IqMsZoN, n=l,...,iVfe. (47) 

M2 

In other words, instead of choosing Mi rows of that span the Mi-dimensional subspace along which 
the energy of signal is maximum as in Section 1111- A I we choose Mi identical rows for each as the 
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1 -dimensional subspace along which the energy of signal is maximum, and then decrease the noise power 
(and therefore increase the effective SNR) by averaging over all Mi measurements. Similarly, for $o^ 
instead of choosing its M 2 rows that span the M 2 -dimensional subspace along which the energy of signal 
is minimum as in Section lTlT-AI we choose M 2 identical rows for $0 each as the 1-dimensional subspace 
along which the energy of signal is minimum, and thus decrease the noise energy (and therefore increase 
the effective SNR) by averaging over all M 2 measurements. 

The following theorem introduces the test for performing ([Hi based on measurement design in (l40l) 
and (|4T1) . 


Theorem 4. Consider the following test for compressive detection of subspace signal Hx 


r = 


Nt 

Ml \zs[n]\^ 
n=l 

Wt 

M2 E koNP 

n=l 


^ 7 


(48) 


The probability of false alarm PpA and the probability of detection Po are then expressed as 

PfA = Q.F{Nb,Nb)^'^)^ 

and 


Pd - QF(N,,Nt)[-^ 


+ M2crlp% 


7 


(49) 


(50) 


+ Mialpl 

Proof: To derive the probability of false alarm PpA and the probability of detection Pd, we first 
notice that if the Gaussian noise is i.i.d, then it is easy to observe that the distribution of Zo[n] is 


ZoU ~ 


A^(0, 


MM 2 


under TLq, 


(51) 


■^(0’ + J 4 ^xP%) under Pi. 


Similarly, the distribution of Zs[n] is 


A /(0 


’ MMi. 


under TLq, 


zAn\ 


(52) 


A" + under Pi, 


Then by following steps similar to those in Theorem |2l (|4^ and dSOl) are concluded. ■ 

The Fully-Correlated Compressive Detector of (l40l) and (Hdl) whose probabilities of false alarm and 
detection were derived in Theorem |4] is optimum in terms of providing the largest difference between 
the expected energy of its two set of samplers. However as discussed earlier, a major measure of interest 
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in statistical hypothesis testing is the power of the test. The following corollary discusses the optimality 
of the proposed Fully-Correlated Compressive design and detector in terms of the power of the test. 


Corollary 2. Among all fully-correlated measurement designs of the form 


and 






1 


y/M 


1 


LMi,1 




•'2 ) 


(53) 


(54) 


VaI 

with ||ti ||2 = ||t 2||2 = 1 <^nd t^t 2 = 0, the one proposed in and d4il) delivers the highest Po for a 
fixed PfA for performing the test given in 


Proof: By taking steps similar to the proof of Theorem H] the probability of false alarm can be 
derived as PfA = QF{Nb,Ni){'y)- Therefore, for a given PfA, Mi, M 2 , and Ni,, we can write 

'(Tg + M2CJ^t^HH'^t2^_i 


Pd — QF{Nb,Nb) 


(Jn 




(Pfa)^ ■ 


(55) 


From t l55D . it is clear that the most powerful design is the one that minimizes which in 

turn is the one that minimizes HH^t 2 and maximizes Then by employing the Rayleigh-Ritz 

theorem / l^ Chapter 8] it can be easily proven that the optimum design is ti = ui and t 2 = u^v- ■ 


Again, similar to the Maximally-Uncorrelated case, Theorem |3] and Corollary |2] result in identical 
solution although they are trying to optimize different cost functions, which is because the proposed 
detector is in fact an F-test detector working based on the difference between the variances of the two 
sets of measurements. 

The performance of Fully Correlated Detector has a clear connection to the way we choose the number 
of samplers. The following Corollary discusses this. 


Corollary 3. For a given number of samplers M and probability of false alarm PfA, the performance 
of the detector Po is maximized when M 2 = 1 and Mi = M — 1. 

Proof: First, we notice from (|4^ that the threshold level 7 is independent of the choice of Mi and 
M 2 . Then from (l50l) it is clear that Po achieves its maximum value when M 2 reaches its minimum value 
and Ml reaches its maximum value. ■ 

The above Corollary asserts that in the case of Fully-Correlated Detector, for a fixed hardware budget, 
in terms of the total number of samplers M = Mid-M 2 , the best performance, in terms of Po, is obtained 
if we employ just one sampler for $0 and spend the rest of the hardware budget on Besides, by 
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Fig. 3. Probability of Detection as a function of Mi for the two detectors introduced in Sections IlII-AI and IIII-BI The total 
number of samplers is M = 10, the Signal-to-Noise Ratio is ^ = OdB and the largest and smallest eigenvalues of HH^ are 
Pi = 10 and p%=0, respectively. The Probability of False Alarm is Pp 4 = 0.01 and the number of temporal measurements is 
Nb = 5 . 


putting the optimum values M 2 = 1 and Mi = M — 1 in (l50b it is clear that the performance enhances 
as the total number of measurement devices M grows. Figure [3] demonstrates an example of how the 
prohahility of detection for Fully-Correlated detector changes with the number of maximum energy 
samplers Mi. The performance curve of the Maximally-Uncorrelated case is also shown for comparison. 
As it can be observed, while the performance of Fully-Correlated detector clearly enhances with increasing 
Ml, the performance of Maximally-Uncorrelated Detector does not show such a monotonic relationship. 
The best performance of the Maximally-Uncorrelated Detector here is achieved for Mi = M 2 = 5. 
The probabilities of detection of the two detectors here have been calculated from (l50l) and (fTVl) . The 
parameters of the Maximally-Uncorrelated detector has been chosen such that Pd = PD,ib = PD,ub- 
Remark 4: Comparing (1^ and (l43l) to (fTTl) and (fT^ it may seem that the performance of the Fully- 
Correlated detector is always superior to the Maximally-Uncorrelated detector as the optimizations in 
(HU) and (l43l) are unconstrained compared to the optimizations in (fTTl) and (fT^ which are constrained. 
However, we should notice that this is not always the case as Maximally-Uncorrelated detector provides 
a higher diversity in its collected samples, which is reflected by the higher degrees of freedom in (fTSl) . 
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([TT]) and (fT^ compared to (|4^ and (l50l) . For example consider the case when x[n] is present hut lies 
orthogonal to the strongest right singular vector of H. Then for the Fully-Correlated samplers the energy 
of samples collected hy all of the Mi samplers includes only the energy of noise while for the 
Maximally-Uncorrelated Compressive Detector only the output energy of the first max-energy sampler 
is equivalent to the noise energy and the output of the rest of Mi — 1 samplers include both the signal 
and the noise energy. This can also he seen from the example of Figure |3] where the best performance 
of the Fully-Correlated detector occurring at Mi = M — 1 = 9 is inferior to the best performance of 
Maximally-Uncorrelated detector occurring at Mi = M/2 = 5. □ 


C. Special Case of Known Noise Variance 

As discussed earlier in this Section, the idea of having two sets of measurement devices and is 
to have two sets of measurements with maximum difference in their variances when the signal truly exists 
and equivalent variances when the signal does not exist and then making decision about the existence of 
signal based on this difference. However, if the variance of noise is known, we can simplify the proposed 
design by discarding as we can simply compare the variance of the output of to the known noise 
variance. In other words, we can employ only Mi measurement devices determined by the rows of 
and then compare the power of measurements taken by them to the power of noise: if they are almost 
the same then we decide that no signal exists, if not then we decide that signal exists. 

Similar to the two detectors proposed in sections ITTT-AI and ITTT-Rl we can have two detectors when 
the noise variance is known. The following two propositions express these two detectors and their 
performances. We remark that since is discarded when noise variance is known, we will have M 2 = 0 
and therefore the whole number of measurement devices is M = Mi. 


Proposition 1. If we design the measurement matrix according to ((91), then in the case of known ctq, the 
test in m can be replaced by 




r = 


n=l 


all Ml 


Hi 

^ 7 - 

-Ho 


The Probability of False Alarm, PfA, for this test is computed as 

PpA = Qx^(MiNt)(7), 

and the Probability of Detection, Pd, is bounded as 

Pd, lb < Pd < PD,ub 


(56) 


(57) 


(58) 
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where PD,ib <ind PD,ub d^e, respectively, the lower bound and the upper bound for the probability of 
detection which will take the following formulas 


PD,lb = Qx'^{MiNt){Vlb'y)^ 

(59) 

PD,ub = Qx^{MiNt){'nub'y)> 

(60) 


where , Puh — and Q^ 2 (d)(x) is the tail probability of the chi-squared distribution 

with parameter (degree of freedom) d at point x. 

Proposition 2. If we design the measurement matrix according to d40l) . then in the case of known Uq, 
the test in d4§l) can be replaced by 


Nt 

E 

/-r' 71 —1 

(Xq/mI ^ 

The Probability of False Alarm, PfA, for this test is 

(61) 

PpA = 

(62) 

and 

(63) 


D. Compressive Detection in Presence of Known Subspace Interference 

Extending the above results to the subspace interference case is straightforward and is similar to the 
same scenario for classical uncompressive detection of deterministic subspace signals in presence of 
interference in iHTI . Suppose that the subspace signal model in ([Til is replaced by 

s[n] = Hx[n] + Gt[n], (64) 

where G is an x it'' matrix and the second term on the right-hand side of (l64l) represents the subspace 
interference. If s[n] were available, the interference could be cancelled by left multiplying (l64l ) by matrix 
= I - G(G^G)-^G^ to obtain 

P^s[n] =P^Hx[n]. (65) 

Comparing (1651 ) to ([T]), it can be easily deduced that in this scenario to design the measurement vectors, 
i.e., rows of the left-singular vectors of matrix H in the interference-free scenario in (l9l) and (ITOl) 
must be replaced by left-singular vectors of matrix PqH right-multiplied by matrix Pq. Otherwise, all 
the detector developments are similar to the previous subsections. 
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IV. Compressive Detector with Non-Ideal Measurements 

In this section we investigate the case when due to the lack of complete knowledge about the suhspace 
or due to the limitations in designing hardware with infinite precision, the measurement vectors do not 
exactly match the singular vectors of matrix H. We consider two models for descrihing this uncertainty. 
The first model, defines the finite precision using additional observation noise. The second model, on 
the other hand, chooses the measurement vector from a set of realizable measurement vectors (hardware) 
which are approximations (or quantized versions) of infinite-precision real-valued singular vectors. 


A. Compressive Detection with Imprecise Measurements 

In this section we analyze the detector when the measurement device is not precise. This imprecision 


in measurements can be modeled as 


i[n] = z[n]-h (5 ^/^e[n], 


( 66 ) 


where z[n] denotes the imprecise compressive measurements, z[n] denotes the basic measurements as 
introduced in Section HIl e[n] ~ AA(0 m,1)Im)> and 6 represents the precision of measurements. Thus the 
imprecision is here modeled through additional observation or measurement noise, as in ll42l . 

Now, if we replace Zs[n] and Zo[n] in the test statistic (fT4l) of Maximally-Uncorrelated detector by 


ZoTi = zn 


l:Mi 


and Zo[n] = z[n] 


Mi+1:M 


respectively, then it can be easily verified that the probability 


of false alarm will be the same as in (fTSl) . but the lower and upper bounds of probability of detection in 
([TT]) and (dUl will be replaced by 


(^o + (^Ip%-M2+i + ^ 


PD.m - ^ ^ 


7 ), 


and 


PD,ut, - ^2 + ^ 2^2 + 7 ), 


(67) 


( 68 ) 


respectively. Comparing (l67]) and (IMI) with (fTTI) and (fT^ . it is clear that the performance degrades as 
the precision of measurement devices 5 decreases. This can be however compensated by improving the 
precision through employing several, say L, identical measurement devices instead of each device and 
then averaging over all measurements. In other words, instead of M measurement devices, we employ 
ML to improve the performance. We call L the hardware budget factor as it indicates the number of 
measurement devices that we employ. 
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In this case, the prohahility of false alarm will he still expressed as in (fTSl) while the lower and upper 
hounds for prohahility of detection are modified as 


^ Mb 


Pnib = 


MiNi,M2Nb 


)( 


-1 


and 


PD,ub = 


2Ni)i 


^0 + + Mb 

+ L^^Ipn + Mb-^ 




—7), 


(69) 


(70) 


(MiNb,M2l^b)y ^2 ^ 2.a2p2 + Mb- 
Comparing (l69l ) and (iTOb with (fTTl) and (fT^ . the hardware budget factor L for achieving a performance 


similar or better than the infinite precision case of 5 ^ = 0 is obtained as 

M 


L > 1 + 


bal' 


(71) 


B. Compressive Detection Using a Fixed Set of Measurement Vectors 

The second model in this section for modeling non-ideal measurement devices is formulated as 
choosing the measurement vectors from a fixed set of unit-norm orthogonal vectors whose members 
are not necessarily singular vectors of the system matrix H. These vectors can be hardware-realizable 
approximation of the singular vectors, e.g. quantized versions of them, or simple down-samplers. Assume 
that M = Ml -|- M 2 measurement vectors should be chosen out oi R > M uncorrelated columns of 
N X R matrix ^ = [r/ji, '02, ■ ■ ■, iI’r] (notice that = I/j). The rows of and $0 then have to 

be chosen from the columns of 

The optimum is the one that maximizes the energy of its output. If we want to follow the Maximally- 
Uncorrelated design of Section ITll-AI then the rows of are the Mi columns of ^ for which ||$sH||i? 
is maximized, i.e., 

= arg max trace(HH^^_A4’®'Xr), (72) 

MC'Jl:\M\=Mi 

^ M 

where TZ = {1,2,... , R} and is the sub-matrix of ^ consisting of columns determined by index set 
M.. Similarly, the rows of $0 can be selected as the M 2 columns of for which ||$sH||i7’ is minimized, 
i.e., 

$0 = arg min trace(HH^’J';\ 4 ^Ti) (73) 

^ Ad 

On the other hand, if we want to follow the Fully-Correlated design of Section Illl-RI then we should 
pick the two columns of say 0^ and 0^, that, respectively, maximize and minimize the signal energy 
it^ 0^HH^0 and then design and respectively, as 
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and 




1 

y/M 


1mi,i <8) 


(74) 


* A 

— 


/m 


1 m 2 ,i ® 'ipj 


(75) 


To see the performance of the detectors in this case, we first provide the detector’s formulation when 
the noise variance is known. It is easy to verify that when the elements of A4 are orthonormal, PpA is 
computed as in (l57]) and (l62l) for the Maximally-Uncorrelated and Fully-Correlated detectors, respectively. 
However, Pd of Fully-Correlated detector as well as the hounds of Pd of Maximally-Uncorrelated detector 
will differ from what we have in (l59l ). (I60l ) and (|6^ . The reason for this is that the rows of are not 
any longer the singular vectors of H and therefore the set of singular values of $<jH is not a subset of 
the set of singular values of H. In other words, (|5^ . (l60l) and (l6^ should he respectively replaced hy 
the following equations 


and 


PD,lh = QxHMiNt)iVlb'l), 
PdM = Qx^{MiNb){Vub7), 


(76) 


(77) 


and 


Pd = 




CTn 


+ Mialpl 


7 


(78) 


where fjif, = ^ Vub — and pi and pMi are the largest and smallest singular values of 

$sH, respectively. Since the eigenvalues of are smaller than the corresponding eigenvalues 

of HH^ |[^ Theorem 6.76], the performance is lower than those indicated in (l59l ). (I60l) and (|6^ . 

When the variance of noise is unknown, on the other hand, we cannot write the prohahility of detection 
of the detector based on P-distribution as the outputs of and measurement devices are not 
independent of each other, which is in turn because the rows of and $o are not any longer the singular 
vectors of H. However, we can still design a constant false alarm rate detector since in the absence of 
signal term Hx the output of these two measurement devices are still independent as = Omi.Ms- 


V. Compressive Detection on a Finite Union of Subspaces 

In this section, we generalize the results of previous sections to the case when we know that the H G S 
where 

S = {Hi,H2 ,...,Hq}, (79) 
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where Hg, g = isa full-column rank N x Ki matrix and the prohahility that the signal lies in 

suhspace Hg is P(H = Hg) = vTg, q = 1,... ,Q. The goal is then to accomplish the test of (lUl. 

Similar to the case of Q = 1, the measurement matrix is tailored to performing the F-test. In other 
words, it must consist of two parts: one lying in the strongest signal suhspace, which attains the maximum 
energy, and the other in the weakest suhspace, which attains the minimum energy. But, since the true q 

is not known a priori, we should design it such that the expected value of the energy over the set S gets 

its maximum value for the first part and its minimum value for the second part. 

To he more precise, if we adopt the Maximally-Uncorrelated strategy, the first and the second sets of 
measurement devices will he chosen as the solutions of the optimization problems 

argmaxEH,x,w(zrzs), s.t. ='^M^, (80) 

and 

argminEH,x,w(zoZo), s.t. = Im^, ( 81 ) 

respectively. On the other hand if we want to design it based on Fully-Correlated strategy we have to 
relax the constraints = Imi and M^o^o — in (1^ and (IMl) . This will guarantee that the 

energy of the signals collected at the first and second channels attain, respectively, the maximum and 
minimum possible values which can then be exploited for optimum compressive signal detection. 


Theorem 5. The first and second channels of measurement matrix which satisfy dSOD and d^il) are 
respectively the strongest Mi eigenvectors and the weakest M 2 eigenvectors of the following matrix 

Q 

Hcq = ^TTgHgH^, (82) 

q=l 


Proof: It is easy to see that 

iii^H,x,w(z5 Z^) 


EH(c72trace(HH^$^$,) + ||) 

Q 2 

^7rgO-2trace(HgH^$J$J + 

q=l 

2 

fj^trace(Heq$f$s) + 


(83) 


Since the second term in the right-hand side of (IMI) is independent of the choice of the measurement 
matrix, the measurement sub-matrix that maximizes (l8^ is the one maximizing the first term, which 
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contains the strongest Mi eigenvectors of Hgq as its rows |[34l P 11.4.5]. The proof of the weakest M 2 
eigenvectors of Hgq being the solution of (IMl) is similar. ■ 

Connection to the compressive detection of sparse signals: The ordinary compressive detection of 
random sparse signals can he specified as a particular case of the setting in this section. To see this, 
assume that we have a iT-sparse signal in A^-dimensional basis as 

s = fix, (84) 

where s E and x E are both vectors of size N and only K N entries of x are nonzero. We 
assume that the K nonzero entries of x are i.i.d with Gaussian distribution. Denoting the true unknown 
support of X by Z, the signal model in (IMl) is equivalent to 


s = fijx, 


(85) 


where x E and fix E contains the columns of fi corresponding to index set Z. In the absence 

of any further structure or knowledge of Z, any K entries of x can be its true support set equally likely 
with probability 1 / . In this case, the members of the union set S in (1791 ) are the Q = size-K 

combinations of the N columns of fi. Furthermore, we will have vri = ... = ttq = ^. 

If the basis fi is orthonormal, then it is easy to verify from (l82l) that 


^eq — 


( 86 ) 


Since all the singular values of Hgq in 


are equal, then Hgq has no structure and therefore the 


n 


but we can 


methods proposed in Sections 1111- A I and Illl-RI are not suitable for detection of this signaf] 
use the method of Section ITll-C I provided that the noise variance is known. However, if there is a known 
structure in the sparse random signal model of ([84] ). for instance if the probability of different iT-size 
combinations are not equal or if fi is not an orthonormal basis, in the sense that not all of the singular 
values of Hgq are identical, then the methods proposed in Sections llll- A I and Illl-R I can be employed for 
detection of signal from compressive measurements. 


VI. Simulation Results 

In this section we study the performance of the proposed compressive detection methods via computer 
simulation studies. 

'Notice that in the absence of noise variance knowledge, none of the methods in the literature can detect the signal. 
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The first experiment is to compare the performance of the two methods proposed in Sections ITTT-AI 
and ITTT-RI which are able to detect the signal under unknown variance condition. In this experiment we 
choose N = 1000 and K = 10. The elements of matrix H are generated independently from Af{0,1/y/K) 
such that the expected value of the squared second norm of its rows is one. Signal x and noise w are 
both drawn from zero-mean Gaussian distribution with variances cj^ and cjq, respectively. We define 
the SNR as lOlog^o ^ (in dB). The number of measurement devices is M G {2,4,8} and we choose 
Afi = M 2 = M/2. The results are obtained from 10000 random trials to gather sufficient statistics for 
reliable evaluation of the performance. 

Figure |4] illustrates the performance in terms of probability of detection as a function of The curves 
have been depicted for two cases of Nf, = 5 and = 50. As it can be observed, while Maximally- 
Uncorrelated compressive detector provides better performance for smaller Nf, at high SNRs, it becomes 
inferior to Fully-Correlated detector for higher A^f,’s and at lower SNRs. This is also clear from Figure |5] 
where the performance in terms of Receiver Operating Characteristics (ROC) has been depicted. As we 
can see, again Maximally-Uncorrelated compressive detector outperforms Fully-Correlated method for 
high SNRs and small A^b’s, while Fully-Correlated method works better in small SNRs and for higher 
number of measurement devices. 

The second experiment is to study the performance loss caused by the lack of knowledge about the 
noise variance. Figures |6] and |7] show the performance of the proposed Maximally-Uncorrelated detection 
method for two cases: when the noise variance is known and when it is unknown. As we expect, the lack 
of knowledge of noise variance results in considerable performance degradation. Flowever, reasonable 
performance is still obtained. 

In the third experiment we study the effect of imprecise measurements and compensating it by deploying 
multiple measurement devices. The imprecision is modeled as described in Section IIV-AI Figure [8] 
illustrates the ROC curves for precise (red curves) and imprecise (blue curves) measurements when 
(7^ = 1, and (7g = 100. The precision of all imprecise devices (blue curves) is chosen as = 100. As 
expected and is clear from the comparison of the red curve with blue curve with L = 1, in equivalent 
number of measurement devices, the imprecision in measurements results in a considerable loss in 
performance. However, this can be compensated by employing higher number of measurement devices. 
For instance, when L = 3, the performance becomes equivalent to the precise case. This can be also 
verified by (17X1) in which the lower bound for performance improvement is L = 3. When we go beyond 
this bound, e.g. L = 5, the performance improves compared to the precise case. 

Figure 0 shows how the probability of detection changes with the noise power for various precisions 
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Fig. 4. Probability of Detection as a function of SNR for Maximally-Uncorrelated compressive detector (Blue) versus Fully- 
Correlated method (Red) when PpA = 0 .05. Solid curves represent the case when Ni, = 50, and dashed curves represent 
the case when Nb = 5. The number of measurement devices is M = 2 (square), M = 4 (circle), and M = 8 (diamond). 
Furthermore we choose Mi = M 2 = M/2. 


and hardware budgets of L G {1,3}. The black curve depicts the performance with a single precise 
device. Here the performance enhancement with deploying L = 3 times devices has been depicted. As it 
can be seen for a certain 6, the performance of the case of L = 3 is better than a single precise device 
as long as ctq < On the other hand, when the variance of noise is higher than this limit, even 

employing L times devices is not enough for compensating the imprecision. 

In the fourth experiment we study the performance of the proposed compressive detection method 
when the signal comes from a union of some subspaces with certain probabilities. For this experiment, the 
ambient dimension is = 20 and the signal dimension is K = 3. In each random trial we first generate an 
N X N matrix with entries from a normal distribution with mean zero and variance 1/K. The elements 
of union set S in (1791) are then Q different iT-size combinations of columns of the generated matrix 
while the signal is generated from one of them which is chosen uniformly randomly from all Q possible 
combinations. This implies that tti = ... = ttq = 1/Q. The experiment is done for Q G {1,10, 50}. We 
remark that the case of Q = 1 coincides with the basic case described in Section [III] as in this case there 
is no uncertainty over the subspace on which the signal lies. 
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Fig. 5. ROC for Maximally-Uncorrelated compressive detector (Blue) versus Fully-Correlated method (Red). Solid curves 
represent the case when Ni, = 50 and o^./oq = —30 dB, and dashed curves represent the case when Ni, = 5 and (y%/o% — —15 
dB. The number of measurement devices is M = 2 (square), M = 4 (circle), and M = 8 (diamond). Furthermore we choose 
Ml = Ms = M/2. 


Figure [TO] illustrates the performance in terms of probability of detection for various values of 
The probability of false alarm here is set to PpA = 0.05 and the number of measurements is Nh = 50. 
The results are shown for M G {2,4} with Mi = M 2 = M/2. Figure [TT] also shows the ROC curves 
when = —lOdS. As it can be seen, as the cardinality of the subspace union set S increases, 

the performance of the algorithm degrades. This is in turn because as the number of possible subspaces 
increases, the rows of measurement matrix $ are less matched to the singular vectors of the true subspace 
matrix. 


VII. Conclusion 

The problem of random signal detection from compressive measurements when the signal is lying (or 
leans toward) a low dimensional subspace was studied. Having the knowledge of the subspace structure, 
we proposed two measurement designs and formulated the hypothesis test and its performance metrics 
for each of them in the case that the noise variance is unknown. We also showed how the design can be 
simplified in the case that the noise variance is known. We analyzed the effects of imprecise measurements 
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Fig. 6. Probability of Detection as a function of SNR for Maximally-Uncorrelated compressive detector under unknown variance 
(blue) versus known variance (red) when PpA = 0 .05. The number of measurements over time is Ni, = 5. The number of rows 
of is Ml — 1 (square). Mi = 2 (circle), and Mi = 4 (diamond). When the variance is known we have M = Mi and when 
it is unknown we choose M 2 = Mi and therefore we have M = 2Mi. 


and showed how it can he compensated hy deploying more identical measurement devices. The problem 
was also generalized for the case when signal belongs to a union of finite number of subspaces with 
known probabilities. 
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